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AE8/AP8  Implementations  in  AE9/AP9, 
IRBEM,  and  SPENVIS 


The  CmdlineAe9Ap9,  SPENVIS  and  IRBEM  implementations  of  the  AE8/AP8  and 
SHIELDOSE2  models  were  used  with  a  GTO  satellite  ephemeris  for  10  full  orbits,  for  the 
calculation  of  electron  and  proton  integral  and  differential  flux  values,  and  their  associated  dose 
values  for  the  three  defined  geometries.  The  flux  values  and  dose  results  from  these  three 
implementations  closely  match,  with  the  recent  revisions  of  the  CmdlineAe9Ap9  (included  in 
Vl.l)  and  IRBEM  library  (to  be  included  in  a  2013  release).  Some  small  differences  are 
observed,  but  can  be  attributed  to  differences  in  the  magnetic  field  models  and  AE8/AP8  model 
settings,  and  for  dose  results  to  differences  in  the  SHIELDOSE2  processing  settings. 

1.  Introduction 

This  document  provides  a  comparison  of  results  from  the  AE8/AP8  model  implemented  within 
the  AE9/AP9  model  application  (CmdLineAe9Ap9),  IRBEM  (previously  called  ONERA)  and 
SPENVIS.  Results  for  time-average  fluxes  and  total  dose  calculations  via  SHIELDOSE2  are 
provided.  The  ephemeris  for  10  full  orbits  of  a  Geosynchronous  Transfer  Orbit  (GTO)  is  used 
for  this  comparison,  as  it  samples  a  wide  range  of  locations  within  the  radiation  belts.  The  new 
AE9/AP9  model  application  includes  the  ability  to  query  the  legacy  AE8  and  AP8  models.  The 
AE8  and  AP8  flux  and  dose  results  from  this  show  some  differences  when  compared  to  those 
from  the  IRBEM  and  SPENVIS  implementations.  For  reference,  the  AE9/AP9/SPM  model  is 
described  by  Ginet  et  al.  [2013]  and  Roth  et  al.  [2014],  AP-8  by  Sawyer  and  Vette  [1976],  AE-8 
by  Vette  [1991],  and  SHIELDOSE2  by  Seltzer  [1994]. 

2.  Generation  of  Model  Results 

The  Ae9Ap9  package  GUI  program  was  used  to  generate  ephemeris  information,  and  perform 
AE8/AP8  and  SFIIELDOSE2  model  runs  via  the  CmdLineAe9Ap9  application.  The  AE8/AP8 
model  is  run  for  solar  maximum,  with  SAA  translation.  SFIIELDOSE2  model  parameters  are 
listed  below.  An  external  script  processed  the  time-tagged  results  to  produce  the  sets  of  average 
integral  and  differential  flux  values  as  a  function  of  energy,  for  comparison  to  the  SPENVIS 
output  products. 

The  information  in  the  CmdLineAe9Ap9-generated  ephemeris  file  was  processed  by  an  external 
script  to  produce  a  file  in  the  fonnat  expected  by  the  SPENVIS  web-based  system;  most 
importantly,  this  script  converted  the  original  Modified  Julian  Date  time  values  to  the  alternate 
definition  used  by  SPENVIS.  This  ephemeris  file  was  processed  by  the  SPENVIS  system  to 
generate  the  AE8/AP8  [solar  max]  integral  and  differential  flux  values,  and  associated 
SHIELDOSE2  model  dose  results  (using  the  same  model  parameters  listed  below). 

The  generated  ephemeris  file  was  also  used  as  input  to  a  program  interfacing  with  the  IRBEM 
library  implementations  of  the  AE8/AP8  models  and  the  SFIIELDOSE2  model.  External  scripts 
were  used  to  produce  sets  of  the  average  flux  values  as  a  function  of  energy. 
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For  a  secondary  comparison,  the  CmdLineAe9Ap9  application  was  also  used  to  produce  an 
equivalent  set  of  flux  and  dose  output  files  from  the  AE9  and  AP9  models,  with  these  results 
included  in  the  comparison  plots  here  for  reference.  Results  from  AE9  and  AP9  are  generally 
somewhat  different  from  results  from  the  various  AE8  and  AP8  implementations.  This  is  to  be 
expected,  given  that  the  AE9/AP9  models  use  new  datasets  covering  a  longer  time  period  than 
prior  models,  and  are  processed  in  a  novel  architecture  providing  additional  capabilities.  The 
AE9  and  AP9  results  shown  are  for  a  specific  time  period  and  type  of  orbit,  and  may  have  a  very 
different  character  in  relation  to  the  AE8  and  AP8  results  for  other  time  periods  and/or  orbit 
types.  These  AE9/AP9  results  are  from  “Mean”  mode  model  runs  only.  The  results  from 
“Perturbed  Mean”  or  “Monte-Carlo”  mode  model  runs  will  also  differ,  as  these  modes 
incorporate  the  uncertainties  due  to  measurement  errors  and  (for  Monte-Carlo  only)  estimates  of 
the  dynamic  variations  due  to  space  weather  processes.  For  more  discussion  of  AE9/AP9 
datasets  and  architecture,  and  the  resulting  differences  from  AE8/AP8,  see  Ginet  et  al.  [2013]. 


3.  Model  Run  Parameters  Used  in  This  Comparison 


Orbital  Ephemeris: 

•  Geosynchronous  Transfer  Orbit  (GTO):  180x35870  km,  inclination=6°,  period=37922 
seconds. 

Eccentricit  0.7314  4,  MeanMotion  2.27  3  0  or  its  da  ,  and  0.0 

•  Duration/Rate  10  or  its,  at  t  37  sec,  100  points  per  or  it. 

•  Time  period:  01  Jan  2013  0000UT  -  05  Jan  2013  0921  UT  (4.386  days)  [arbitrarily  chosen] 
(Note:  this  time  period  applies  to  orbit  generation  only;  AE8/AP8  utilizes  geomagnetic  field 
models  from  other  epochs  as  specified  in  the  table  below.) 


SHIELDOSE2  model: 

•  Materials:  aluminum  shielding;  silicon  detector 

•  Dose  depths:  0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.8,  1,  1.5,  2,  2.5,  3,  4,  5,  6,  7,  8,  9,  10,  12,  14,  16,  18, 
20  mm 

•  Geometries:  Spherical,  Finite  Slab,  Semi-Infinite  Slab 

•  Nuclear  attenuation:  none 


4.  Discussion  of  Differences  between  Model  Implementations 


Differences  between  the  results  from  CmdLineAe9Ap9,  IRBEM  and  SPENVIS  are  generally 
less  than  5-40%.  These  differences  are  due  to  several  issues: 

•  The  magnetic  field  models  for  the  implementations  are  different,  as  summarized  in  Table  1. 
Neither  IRBEM  nor  SPENVIS  use  an  external  field  model.  Most  of  the  differences  in  fluxes 
are  explained  by  the  differences  in  the  default  internal  field  model.  With  appropriate 
advanced  settings  in  AE8/AP8,  the  IGRF  field  model  closely  matches  the  Jensen  and  Cain 
field  model  and  consequently  output  fluxes  are  very  similar.  The  lack  of  an  external  field  in 
SPENVIS  and  IRBEM  can  result  in  significant  differences  at  high  L-values. 
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Table  1 


Application 

AE8MIN 

AE8MAX 

AP8MIN 

AP8MAX 

SPENVIS 

Jensen  and  Cain 

Jensen  and  Cain 

Jensen  and  Cain 

GSFC  12/66, 
updated  to  1970 

IRBEM 

Jensen  and  Cain 

Jensen  and  Cain 

Jensen  and  Cain 

GSFC  12/66, 
updated  to  1970 

CmdLine 

Ae9Ap9 

IGRF  @ 
date/1962 

Ext:  Olsen-Pfitzer 

IGRF  @ 
date/1962 

Ext:  Olsen-Pfitzer 

IGRF  @ 
date/1962 

Ext:  Olsen-Pfitzer 

IGRF  @ 
date/1970 

Ext:  Olsen-Pfitzer 

w/  SAA  Translation 

Rotated  0.3°/year  from 

1962  to  run  year 

Rotated  0.3°/year  from 

1962  to  run  year 

Rotated  0.3°/year  from 

1962  to  run  year 

Rotated  0.3°/year  from 

1970  to  run  year 

AE9 

AP9 

CmdLine 

Ae9Ap9 

IGRF  @  appropriate  date 

Ext:  Olsen-Pfitzer 

IGRF  @  appropriate  date 

Ext:  Olsen-Pfitzer 

References :  Jensen  and  Cain  [1962];  GSFC  12/66,  updated  to  1970  [Cain  et  al.,  1967]; 

IGRF  (no  extrapolation  of  field  beyond  01  Jan  2015)  [IAGA,  2010];  Olsen-Pfitzer  [1977] 


•  The  interpolation  methods  applied  to  the  AE8/AP8  flux  maps  are  different.  For  interpolation 
in  spatial  location,  the  implementation  within  CmdLineAe9Ap9  uses  the  original  NASA- 
published  algorithm,  whereas  the  SPENVIS  uses  an  improved  version  for  smoother  maps 
[Heynderickx  et  al.,  1996].  This  can  result  in  differences  between  the  profiles  of  flux  vs. 
time  or  location  between  the  implementations,  with  the  CmdLineAe9Ap9  implementation 
results  being  less  smooth  (however,  overall  fluxes  are  still  comparable).  In  addition,  the 
original  AE8/AP8  models  provided  only  integral  flux  values.  Their  subsequent  conversion  to 
differential  flux  may  be  implementation-specific  as  well. 

•  The  application  defaults  of  the  respective  AE8/AP8  model  implementations  are  different,  so 
proper  comparisons  will  require  adjustment  to  these  settings.  In  CmdLineAe9Ap9,  the 
‘Advanced’  magnetic  field  option  ‘SAA  Translations’  should  be  set  to  ON  to  more  closely 
match  the  magnetic  field  model  used  in  SPENVIS  and  IRBEM,  and  therefore  produce  a 
better  match  of  the  flux  values. 

•  The  SPENVIS  AE8/AP8  model  implementation  uses  fixed  sets  of  30  energy  levels  for 
electrons  and  protons,  while  the  CmdLineAe9Ap9  implementation  uses  fixed  sets  of  2 1 
energy  levels  each.  Some  of  the  differences  observed  in  the  differential  flux  results  may  be 
partly  due  to  the  relative  spacing  of  the  energy  levels  within  their  respective  sets. 

•  The  fonn  of  the  input  ephemeris  time  values  is  different  between  these  three  applications. 
IRBEM  routines  require  a  year,  day  of  year,  and  seconds  of  day.  The  CmdLineAe9Ap9 
application  uses  Modified  Julian  Date  for  input  (and  output)  of  time  values.  SPENVIS  is 
also  documented  to  use  “Modified  Julian  Date”,  but  uses  a  different  definition  for  this  value 
(see  Alerts  to  Users  below). 

5.  Recent  Modifications  of  the  CmdLineAe9Ap9  and  IRBEM 
Implementations 

•  CmdlineAe9Ap9:  The  routines  for  AE8  and  AP8  models  have  been  revised  to  calculate  the 
differential  flux  values  from  the  integral  flux  values  using  a  power-law  method,  instead  of  a 
linear  method.  This  correction  is  implemented  in  AE9/AP9/SPM  V 1 . 1 . 
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•  CmdlineAe9Ap9:  The  SHIELD0SE2  model  database  containing  the  Bremsstrahlung  data 
tables  has  been  updated  to  correct  the  error  present  in  the  original  SHIELDOSE2  publication. 
The  error  is  in  the  ‘elbrbas2.dat’  file:  except  for  A1  detector  targets,  the  finite  and  semi¬ 
infinite  slab  data  tables  for  Bremsstrahlung  have  been  reversed  [Heynderickx,  private 
communication,  May  2013].  This  correction  is  implemented  in  AE9/AP9/SPM  VI .  1 . 

•  IRBEM  library:  The  routine  used  for  the  (internal)  loading  of  the  Bremsstrahlung  data  tables 
has  been  modified  to  correct  the  error  present  in  the  original  SHIELDOSE2  publication 

[ Heynderickx ,  private  communication,  May  2013].  This  correction  is  implemented  in 
IRBEM  library  as  of  July  2013. 

6.  Alerts  to  SPENVIS  Users 

•  SPENVIS  uses  a  non-standard  definition  of  Modified  Julian  Date  (MJD)  for  date  input.  MJD 
is  normally  defined  as  JD-2400000.5,  such  that  0000  GMT  on  17  Nov  1858  is  MJD=0.  In 
SPENVIS,  MJD  is  defined  (for  input)  as  time  from  0000  GMT  on  01  Jan  1950  (see 
http://www.spenvis.oma.be/help/models/sapre.html). 

•  The  implementation  of  SHIELDOSE2  in  the  CmdLineAe9Ap9  application  and  the  IRBEM 
library  follows  the  original  SHIELDOSE2  code  in  that  doses  for  spherical  geometries  are 
calculated  for  fluxes  over  a  hemisphere,  rather  than  from  a  full  sphere.  Therefore,  these  dose 
results  should  be  multiplied  by  2  when  comparing  them  to  the  full-sphere  dose  results  from 
SPENVIS. 


7.  Plots 
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Proton  Fluxes: 


Integral  Fluxes  are  a  near-exact  match  between  the  three  AP8  model  implementations. 

Differential  Fluxes  show  some  minor  differences,  attributable  to  differences  in  the  magnetic  field 
model  and  processing  settings. 


Protons:  AP8  Integral  Flux  Calculations 


Protons:  AP8  Differential  Flux  Calculations 
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Electron  Fluxes: 


Integral  Fluxes  are  a  near-exact  match  between  the  three  AE8  model  implementations. 

Differential  Fluxes  show  some  minor  differences,  attributable  to  differences  in  the  magnetic  field 
model  and  processing  settings. 


Electrons:  Ax8  Integral  Flux  Calculations 


Electrons:  Ax8  Differential  Flux  Calculations 
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Differential  Rux  Ratios  Differential  Rux  Ratios 


Ratios  of  Differential  Fluxes: 


Differential  fluxes,  used  as  input  for  the  dose  calculations,  show  some  minor  differences  relative 
to  SPENVIS,  attributable  to  differences  in  the  magnetic  field  model  and  processing  settings. 


Ratios  ofAx8  Differential  Flux:  Protons 


Ratios  ofAx8  Differential  Flux:  Electrons 
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Proton  Dose  Results:  Finite  Slab 

(for  protons,  the  results  for  Semi-infinite  slab  geometry  are  identical) 

Doses  are  mostly  matching,  with  some  minor  differences  that  are  attributable  to  the  small 
differences  in  the  input  differential  flux  values.  Use  of  the  IRBEM  implementation  of 
SHIELDOSE2  with  the  CmdlineAe9Ap9  application  and  SPENVIS  differential  flux  values 
reproduce  their  respective  dose  results. 


Protons:  Ax8  total  dose  in  Finite  Slab 


Ratios  of  SD2  Proton  Dose  in  Finite  Slab  from  Ax8  Differential  Flux 
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Proton  Dose  Results:  Full  Sphere 

Doses  are  mostly  matching,  with  some  minor  differences  that  are  attributable  to  the  small 
differences  in  the  input  differential  flux  values.  Use  of  the  IRBEM  implementation  of 
SHIELDOSE2  with  the  CmdlineAe9Ap9  application  and  SPENVIS  differential  flux  values 
reproduce  their  respective  dose  results. 


Protons:  Ax8  total  dose  in  Sphere 


Ratios  of  SD2  Proton  Dose  in  Full  Sphere  tom  Ax8  Differential  Rux 
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Electron+  Bremsstrahlung  Dose  Results:  Finite  Slab 

Doses  match  very  closely,  with  some  minor  differences  that  are  attributable  to  the  small 
differences  in  the  input  differential  flux  values.  Use  of  the  IRBEM  implementation  of 
SHIELDOSE2  with  the  CmdlineAe9Ap9  application  and  SPENVIS  differential  flux  values 
reproduce  their  respective  dose  results. 


Electrons  +■  Bremsstrahlung:  Ax8  total  dose  in  Finite  Slab 


Ratios  of  SD2  Electrons+Eremsstrahlung  Dose  in  Finite  Slab  from  Ax8  Differential  Flux 
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Electron+  Bremsstrahlung  Dose  Results:  Semi-Infinite  Slab 

Doses  match  very  closely,  with  some  minor  differences  that  are  attributable  to  the  small 
differences  in  the  input  differential  flux  values.  Use  of  the  IRBEM  implementation  of 
SHIELDOSE2  with  the  CmdlineAe9Ap9  application  and  SPENVIS  differential  flux  values 
reproduce  their  respective  dose  results. 


Electrons  +  Bremsstrahlung:  Ax8  total  dose  in  Semi-Infinite  Slab 


Ratios  of  SD2  Electrons-H3remsstrahlung  Dose  in  Semi-Infinite  Slab  from  Ax8  Differential  Flux 
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Electron+  Bremsstrahlung  Dose  Results:  Full  Sphere 

Doses  match  very  closely,  with  some  minor  differences  that  are  attributable  to  the  small 
differences  in  the  input  differential  flux  values.  Use  of  the  IRBEM  implementation  of 
SHIELDOSE2  with  the  CmdlineAe9Ap9  application  and  SPENVIS  differential  flux  values 
reproduce  their  respective  dose  results. 


Electrons  +■  Bremsstrahlung:  Ax8  total  dose  in  Full  Sphere 


Ratios  of  SD2  Electrons+Bremsstrahlung  Dose  in  Full  Sphere  from  Ax8  Differential  Flux 

2.0 

1.8 

1.6 

1.4 

1.2 
C/3 

O 
TO 

£  1.0 
d) 

CO 

o 
Q 

0.8 

0.6 

0.4 

0.2 

0.0 

0.1  1  10  100 

Shielding  Thickness  (mm) 


IRBEM/SPENVIS  — i — 
CmdLine/SPENVIS  — 
lrbemSD2(SpenvisDF)/SPENVIS  a 
lrbemSD2(CmdLineDF)/CmdLine  v 


Approved  for  public  release;  distribution  is  unlimited. 

12 


8.  References 


Cain,  J.  C.,  S.  J.  Hendricks,  R.  A.  Langel,  and  W.  V.  Hudson  (1967),  A  proposed  model  for  the 
international  geomagnetic  reference  field,  1 965,  J.  Geomag.  Geoelectr.,  19,  p.  335. 

Ginet,  G.  P.,  et  al.  (2013),  AE9/AP9/SPM,  New  models  for  specifying  the  trapped  energetic 
particle  and  space  plasma  environment,  Sp.  Sci.  Rev.,  doi:  10.1 007/s  1 12 14-0 1 3-9964-y. 

Heynderickx,  D.,  J.  Lemaire,  E.  J.  Daly,  and  H.  D.  R.  Evans  (1996),  Calculating  low-altitude 
trapped  particle  fluxes  with  the  NASA  models  AP-8  and  AE-8,  Radiat.  Meas.,  26,  pp.  947-952. 

International  Association  of  Geomagnetism  and  Aeronomy,  Working  Group  V-MOD  (2010), 
International  Geomagnetic  Reference  Field:  The  eleventh  generation,  Geophys.  J.  Int’l., 

183(3),  pp.  1216-1230. 

Jensen,  D.  C.,  and  J.  C.  Cain  (1962),  An  interim  geomagnetic  field,  J.  Geophys.  Res.,  67,  p.  3568. 

Olson,  W.  P.  and  K.  A.  Pfitzer  (1977),  Magnetospheric  Magnetic  Field  Modeling,  Annual 
Scientific  Report,  Air  Force  Office  of  Scientific  Research,  Contract  #F44620-75-C-0033, 
McDonnell  Douglas  Astronautics  Co.,  Huntington  Beach,  CA. 

Roth,  C.  J.,  et  al.  (2014),  AE9/AP9/SPM  Radiation  Enviomment  Model  User’s  Guide, 
AFRL-RV-PS-TR-20 14-00 13,  Atmospheric  and  Environmental  Research,  Inc.,  Lexington,  MA. 


Sawyer,  D.  M.  and  J.  I.  Vette  (1976),  AP-8  Trapped  Proton  Model  Environment  for  Solar 
Maximum  and  Solar  Minimum,  NSSDC/WDC-A-R&S  76-06,  Natl.  Space  Science  Data  Center, 
Greenbelt,  MD. 

Seltzer,  S.  M.  (1994),  Updated  calculations  for  routine  space-shielding  radiation  dose 
estimates,  SHIELDOSE2,  National  Institute  of  Standards  and  Technology  report  NISTIR 
5477. 

Vette,  J.  I.  (1991),  The  AE-8  Trapped  Electron  Model  Environment,  NSSDC/WDC-A-R&S  91- 
24,  NASA  Goddard  Space  Flight  Center,  Greenbelt,  MD. 


Approved  for  public  release;  distribution  is  unlimited. 

13 


DISTRIBUTION  LIST 


DTIC/OCP 

8725  John  J.  Kingman  Rd,  Suite  0944 
Ft  Belvoir,  VA  22060-6218  Icy 

AFRL/RVIL 

Kirtland  AFB,  NM  87117-5776  2  cys 

Official  Record  Copy 

AFRL/RVBXR/Dr.  Michael  Starks  1  cy 


Approved  for  public  release;  distribution  is  unlimited. 

14 


